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1 Introduction 



The structure functions of deep inelastic lepton-proton scattering (DIS) are 
related with an absorptive part of the matrix elements 

i J d 4 xe iqx (p\TJ* m (x)J* m (0)\p). (1) 

Here J^ m {x) is an electromagnetic current and \p) means a proton state. 
In its turn, an expansion of a chronological product of two electromagnetic 
currents near the light-cone in terms of composite operators O a ' m looks like: 

•m 

fib' 

a even m 

+EE ^(^^^^•••^ ro °^,..M m ( ^ 2 )+--- > ( 2 ) 

a even m 

where a = NS, S, g. The dots denote other Lorentz structure, constant term 
as well as gradient terms. The latter give no contribution to DIS structure 
functions. The index a runs all allowed types of the composite operators. The 
expansion (J2J) is a generalization of the operator product expansion (OPE) 
at small distances pQ. It is known as light-cone OPE [2] (see also [3J Hj). The 
quantity Ca}m is called an OPE coefficient function (CF). The composite 
operators in Eq. (J2J) need a regularization. After a renormalization of the 
operators, there arises a dependence of CF's on a renormalization scale \x. 

On the other hand, the DIS structure functions can be represented in a 
factorizable form jH] . For instance, we have 



i 




Here f a {x, M 2 ) is a distribution of a parton of type a inside the nucleon, while 
quantities C a (x,Q 2 /M 2 ) are known as DIS coefficient functions. In Eq. (jHJ), 
M is a factorization scale. It is usually identified with the renormalization 
scale \i. The m-th moment of the DIS coefficient function, C a (m,Q 2 / fi 2 ), is 
identified with the Fourier transform of the corresponding OPE CF's in (0). 
The parton distributions f a {x,ii 2 ) are related with matrix elements of the 
composite operators between one-nucleon states (for details, see 
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The coefficient functions C a (x,Q 2 /M 2 ) were calculated in perturbative 
QCD. For example, contributions of the first and second orders in a s can be 
found in Refs. jH] (where light quarks were considered), and in Refs. [7] (where 
both light and heavy quarks were accounted for). Note that a very definition 
of DIS coefficient functions applies to the diagram (perturbative) technique, 
not to the operator formalism. Moreover, since nucleon wave function is yet 
unknown, one has to deal with diagrams which describe a lepton scattering 
off (nonphysical) quark or gluon off-shell state. 

The light-cone OPE for the scalar theory was studied in Refs. jHj in which 
rules for calculating CF's were presented. Note that the T-product of two 
scalar currents near the light-cone was defined in term of so-called bi-local 
light-ray composite fields |S] . The local light-cone expansion can be obtained 
by performing a Taylor expansion of the non-local one The non-local 
expansion is more general, but we restrict ourselves to considering local OPE. 
In Refs. [10J a problem of finding n-loop contributions to the OPE CF's were 
reduced to evaluating of propagator type (n + l)-loop Feynman diagrams. 

The goal of the present paper is to present a derivation of a closed repre- 
sentation for the OPE CF's in term of composite operator Green functions 
which does not lean on the perturbation theory. This will be done in the 
next Section. In Section |3] we check a validity of our results in a free scalar 
field theory. In Section |U we calculate nonsinglet CF's in perturbative QCD 
in order to demonstrate that our main formulae not only reproduces well- 
known expressions for the quark CF's, but enables us to obtain CF's of the 
gradient operators in the OPE. The finite renormalization of CF's is consid- 
ered in Section 03 In Appendix A a number of useful mathematical formulae 
is collected. In Appendix B we show that our scheme does result in a set of 
homogeneous renormalization group equations for the OPE CF's. 

2 OPE coefficient functions and matrix 
elements of composite operators 

Let us define a quark electromagnetic current 

j; m (x) = V(x)^Q*(x), (4) 
where *ff(x) is a quark field. The electric charge operator in (J3J, 

Q = ^(A 3 + ^=A 8 ), (5) 
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obeys the equations: 



Sp(Q 2 A a ) = ~ + -^=5a8 + 2 5^ . (7) 

Here A a (a = 1, 2, . . . 8) are the Gell-Mann matrices, Sp(X a X b ) = 25 ab , and A 
is the identity matrix. 

The operator product expansion (OPE) for the T-product of two electro- 
magnetic currents looks like (see, for instance, jlj) 



1 00 m -m 

m=0 1=1 



+ ■■-, (8) 



where the dots denote contributions from other Lorentz structures and singlet 
quark and gluon operators. It is clear from Q that this expansion should 
contain nonsinglet (triplet and octet) and singlet composite operators. 

Near the light-cone, a leading contribution comes from twist-2 operators. 
For instance, quark twist-2 (traceless) operator is of the form (operator S 
means a complete symmetrization in Lorentz indices): 

O a Ns,i...»Jx) = i m - 1 Sd, l+1 ...d tlr J(x) 1 , 1 D tl2 ...D, l X^(x) 

+ (terms proportional to (9) 

Here 

D^ = d, + igt a Al (t a = X a /2) (10) 

is a covariant derivative, and A a (x) is a gluon field. 

If the OPE (JHJ) is applied to deep inelastic scattering, only operators of 
the type 

0^ M ,.. M >) = i m - 1 S*(x) 1 , 1 D^...D fMm X a *(x) 

+ (terms proportional to g^^) (11) 

are important, since forward matrix elements of the operators with 1 ^ I ^ 
m—1 are zero. 1 For non- forward Compton scattering, all operators contribute 

1 It is clear from the relation (p | d^O \p') ~ (p — p')^. 
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proportionally to (p — • • • (p ~ v')^ m - The invariant structure which 

survives at p' — > p can be related to the "skew" parton distributions. In our 



notation, O a '"' 



■\a, m,m 

' NS, fli...fl m ~ NS, ft\...ftm' 

In general case, all operators with 1 ^ I ^ m should be take into account, 
since they are mixed under the renormalization, and a renormalized operator 
O^ 1 is defined via unrenormalized operators O™ ' by the relation (we have 



dropped non- relevant indices): 



O 



in, I 



(12) 



where Z is a triangle matrix. In particular, we find that the composite 
operator O m ' 1 is multiplicatively renormalized, while the composite operator 
O m,m , which is relevant to DIS, is not. In perturbative QCD in the first 
order of a s = g 2 /(4n), elements of the matrix Z are given by the following 
expressions A\ (with non-significant finite terms omitted): 



Zw 



a. 



e Vm 2 



Z7T 



e \fi 2 




l-V l+l 



V = I, 

V < I, 
I' > L 



(13) 

where Cp = (N c — 1)/2N C , with iV c being a number of colors. In deriving 
()13|). a dimensional regularization JT] was used, and e = 2 — D/2, where D 
is the number of dimensions. Here and below p, denotes the regularization 
mass, 2 while \x denotes the renormalization scale. 

In what follows, we will be interested in one of nonsinglet quark operators, 



namely, in O 



3, m,l 

NS,Hi...Hr, 



Let us introduce brief notations: 



O™' 1 u (x) 



^3,771,/ 

= U NS,m... 



[*), 



(14) 
(15) 



2 The regularization scale /i arises in dimensional regularization, when one changes an 
integration volume, d A k — ► p, 4 ~ D d D k. 

3 It is obviously that the operator O^}™' (x) can be treated in the same man- 
ner. The singlet case is technically more complicated, since the singlet quark and gluon 
operators have to be considered simultaneously. 



5 



After multiplying both sides of Eq. (jSJ) by the operator 0™f Un (with arbitrary 
n, and k = l,2...n), we get the following relation for T-products of the 
composite operators between the vacuum states: 4 



d'xe^ J d 4 ze^(Tj; m (x)r u m (0) 



m=0 1=1 



x J dW^TC^jO) 0:t^)) + ... . (16) 

Here dots mean other Lorentz structures. It is necessary to stress that prop- 
agator of only singlet quark operator has to appeared in (fTBj) . due to rela- 
tion (J7j). 

As usual, WG ctSSU.1116 t licit Crn I 

(x 2 ) are tempered generalized functions 
(this is explicit in perturbative calculations) so the symbolic relation 



III 



. . . **» = (-2,)- H -j^r (-?T [ w ) (17) 

holds in connection with the Fourier transform in (|T^|) . Then we obtain that 
Vie" 1 J rf 4 ze^(TJ^(x) j; m (0) OJ&Jz)) 

oo m ... .. 

m=o 1=1 v y / 
x J d 4 ze^(TO^j0) O n v t,M) + ■ ■ ■ , (18) 

where C m ,i{.<f) is a Fourier transform of C m j(x 2 ): 5 

Cm,i(q 2 ) = ^ {~q 2 ) m (J^y jd'xe^C^x 2 ). (19) 
Let to be a light-cone 4- vector which is not orthogonal to 4-momentum 

n J = 0, pn 0. (20) 

4 We used the relation T(TOi(.ti) 2 (x 2 )) 3 (x 3 ) = TO^Xi) 2 {x 2 ) 3 {x 3 ). 
5 Some authors include the factor 1/ml in a definition of a composite operator. 
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Throughout the paper, we will work in the limit 

V\ - 0, v\ < 0. (21) 
Let us now convolute our matrix elements with the projector 

(22) 



(pn) n 

In particular, we can define the following invariant structure, 
n Ul ... n v 



= - l - g ^\u,Q 2 ) p 2 ) + ... , (23) 
which depends on variables p 2 , Q 2 = —q 2 , and 



» = - = % (*) 

x Q 2 



The propagator of the composite operator has the following Lorentz struc- 
ture: 



d±z^ z {TO™l^)OltvM) 
= Zp^...p» m Pv i ...Pv n (O m > l O n > k ){p 2 ) 

+ (terms proportional to g^ H p 2 , g^u.p 2 , Q^p 2 )- (25) 
Equation (J25|l means: 

^ 4 ^^(To; ; . , , ra (o)o:i„( 2 )) 



(pn) n 

= 2p^... Pflm (O m > l O n > k )(p 2 ). (26) 

Note that F n ' k (uj,Q 2 } p 2 ) and (O m ' l O n ' h }(p 2 ) are dimensionless. 

Let us note that at p 2 — > propagators of composite operators of higher 
twists are suppressed by powers of p 2 with respect to the propagators of 
twist-2 operators. Thus, our approach enables us to isolate a contribution 
from twist-2 operators. 
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At fixed Q 2 and p 2 , 3-point Green function (T J^ m J^ m O n,k ) has a discon- 
tinuity in the variable (q + p) 2 for (g + p) 2 ^ (that is, for u> ^ 1). By using 
the unsubtracted dispersion relation for F n ' k (uj,Q 2 ,p 2 ), 



F n ' k (oo,Q 2 ,p 2 " 



duo' 



hxxF n ' k {u',q 2 ,p 2y 



it I w — UJ 



oo 

— y 



oo 



m=0 



duo UJ 



/, ,/— m— 1 



disc^'V^V), (27) 



one can derive from Eqs. (fTHjl and (J23*j) . (|231) : 

m 

^^(Q 2 )(O m ''0"' fe )(^ 



«=1 



p 2 ^0 



1 

2vri 



dxx m - 1 disc (p+g)2 F n ' fe (x, Q 2 , p 2 ) 



(28) 



As we will see in the next Sections, both propagator of the composite opera- 
tor and matrix element F n,k need a renormalization already in zero order in 
strong coupling constant [T2j. Thus, we have (p m,l O n ' k ) = (O m,l O n,k )(p 2 / fi 2 ) , 
pn,k _ F n ' k (x,Q 2 /p 2 ,p 2 //i 2 ), and, consequently, C m j = C m j(Q 2 /fi 2 ), where 
p, is a renormalization scale. Remember that all these quantities are dimen- 
sionless. 

The Eq. (I28|) are valid for all integer k ^ 1 (of course, k ^ n, but we can 
choose O n ' fc with any n). Let us put n ^ m and define the matrixes: 



(00) ifc = (TO m > l O n > k )(p 2 /n 



and 



(JJO) fe = (TJ em J em O n > k )(x,Q 2 /p 2 ,p 2 /ii 2 ) 



p 2 ^0 



(29) 



(30) 



where Lorentz invariant part of 3-point Green function is implied. 6 Then 
from (|28|) we obtain m equations for CF's (for k — 1, 2, . . . m): 

d m,i = E h I dxxm ~ X &s C(p+q)2 {JJO) k (00)u\ (31) 

J 1 ** 







6 Thcre is no dependence of IQOjik on indices m, n, and no dependence of (JJO) k on 
index n except for a trivial factor (—1)™. 
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where (OO)^ 1 is the inverse of the matrix (00) ik- 

The formulae (J28|) . (fBTj) is our main theoretical result. The equation ()31|) 
gives an operator definition of the OPE CF's in term of vacuum matrix 
elements of the composite operators. 7 It is important to stress that our 
definition does not lean on a notion of quark distributions. 



3 Coefficient functions in free scalar field 
theory 

As a simple example, let us consider an expansion of T-product of two com- 
posite operators <fi 2 (z) =:<fi(z)<p(z): in free scalar field theory. The scalar 
fields 4>(z) are assumed to be real and massless. Using the Wick theorem, 
one can find 

T0 2 (x) <f> 2 (y) = [D c (x - y)f + 2D c (x - y) : </>(y)(f>(x) : , (32) 

where 

DJz) = — - — (33) 

is a causal propagator of the field 4>(z). By expanding bilocal operator in 
powers of variable (x — y) and putting then y = 0, we get (here and below a 
constant term is omitted): 



■l) m d d d 



T0 2 (x)0 2 (O) = 2D c (x)S" — x^x^.-.x?" 

z — ' ml 

+ a Ml ^ 2 ... ^ m ]0 2 (O). (34) 

As usual, the symbol <9 M ( d^) in (|34|) denotes a derivative which acts on the 
left (right) standing field <fi in the composite operator </> 2 (0). 

Using an explicit symmetry of x^ l x^ 2 . . . of 1 " 1 in indices, one can write 

(-l) m x^ . . . of™ V w V M2 . . . V Mra = i w . . . ^ 

(-1)' I , ) d, l+1 d, l+2 • • • cV 9 W 9 W . . . Ml , (35) 
z=o ^ ' 



7 The electromagnetic current Q is a particular case of a quark composite operator 
with zero anomalous dimension. 
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where d^= + <9 M is a total derivative of (fi 2 (z). Then we obtain from 
Eqs. (JMJ), dnSI): 



(fW 9 * T0 2 (x) 2 (O) 



-g 2 + zO 



m,l 9ml- 



(36) 



m=0 Z=0 

Here (— q 2 + z0) _1 is a Fourier transform of D c (x), and a brief notation, 



ml 



(37) 



is introduced. The OPE CF's in Eq. (J3fij) are the following dimensionless 
quantities: 



n 

^m,rr, 
Cm, i 



1 = 0,1, 



,771 — 1. 



(38) 



Let us demonstrate that our main formula (|28j) gives the same result 
To do this, one has to calculate both the propagator of the composite oper- 
ator, (<fi m i <fin k) i an d 3-point Green function ((fi 2 <fi 2 4>n k) ■ The corresponding 
diagrams are presented in Fig. The vertex which corresponds to the com- 
posite operator 2 k is equal to {pn) n ~ k (kn) k , with k and (k—p) being ingoing 
and outgoing 4-momentum, respectively. The result of our calculations is the 
following: 



(<fiml<filk)(j> 2 /fi 2 



■1) 



71 — 1 



16vr 2 



In 



—p 1 



B(l + l,k+l) 



(39) 



and 



disc 



(p+q) 



p- 1 — >o 



(40) 



The quantity B(x,y) is the beta-function, B(x,y) = B(y,x). 

In order to get a finite result for free scalar propagator (JBTJj) . we had 
to make a renormalization. As is well known, any Green function with an 
insertion of one composite operator is multiplicatively renormalized [T3*l IT1] . 
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The renormalization of Green functions with the insertion of two ( or more ) 
composite operators needs additive counterterms ^2]- The details can be 
found in Appendix B. 

Form formulae (|39jh (|40|) and (j2Hj) (taking into account the factor in front 
of the sum in Eq. (p?B|)). we obtain a set of equations for CF's: 



It is easy to check that C-numbers C m j do obey equations (j4~Tj) for all 
m, k ^ (for this purpose, formula (jA.l|) from Appendix A should be used). 

4 Calculations of coefficient functions in 
perturbative QCD 

In this section we will use our formula (J28)) for calculations of the OPE CF's 
in perturbative QCD. For the composite operator, we will often use a brief 
notation O m ' 1 instead of O 7 ^ 1 „ m (0). Contrary to the matrix element of O m ' 1 
between one-particle (quark) states, {p\O m,l \p}, which defines nonsinglet 
quark distributions in DIS, a propagator (O m,l O n,k ) is divergent already in 
zero order in the strong coupling ot s . The matrix element (J J O n,k ) (Q 2 , x, p 2 ) 
is also divergent, whereas its discontinuity in variable (p + q) 2 is not. The 
2-point Green function (O m,l (x)O n ' h (y)), is renormalized by a subtraction 
of a contact term of the form (Z~ x )\, (Z' 1 )^, f l B k 5(x — y), where Z is the 
renormalization matrix of the composite operators. The details can be found 
in Appendix B. 

We work in the dimensional regularization and use the MS-scheme to 
renormalize ultra-violet divergences. Although all calculations will be done in 
the Feynman gauge, our results are gauge invariant since we sum all diagrams 
in each order of perturbation theory. Remember that in order to find the OPE 
CF's, we have to retain only leading terms in the limit p 2 — > 0. This simplifies 
our calculations significantly. We will restrict ourselves by considering leading 
terms in ln(Q 2 /// 2 ), although our formula (j2Bj) enables one to calculate sub- 
leading terms as well. In other words, along with the limit p 2 — > 0, we are 
interested in large values of Q 2 . 




1 



m + k + 1 



(41) 
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Let us start from the leading (zero) order in a s . 8 The corresponding 
Feynman diagrams are shown in Fig. 0] All diagrams have logarithmic sin- 
gularity at p 2 — > 0. They can be easily calculated with the use of Feynman 
rules for the composite operators presented in Fig. |21 Note that a vertex 
which corresponds to the composite operator is not symmetric with respect 
to ingoing and outgoing momenta. For a particular case with k = n, p = 0, 
we reproduce the well-known expressions (see the first paper in Refs. |15j). 

The result of our calculations in the leading order in the coupling constant 
is the following: 9 

(O m > l O n ' k )W(p 2 /jl 2 ) = (-1)^-1- hi (J*\ B (l + l,k+ 1), (42) 

2it z \ —p z J 

and 

disc (p+(?)2 (JJO n ^°\x,Q 2 /p 2 ) 



(43) 



where e 2 is an electric charge of a quark inside a loop. The terms in square 
brackets in Eq. (|43j) correspond to two bottom diagrams in Fig. 0] with op- 
posite directions of the quark momentum inside the loop. As for two top 
diagrams in Fig. they are identical. So, only one of these diagrams was 
taken into account in Eq. (|42|) . The formula (J42)) gives a finite (singular at 
p 2 — > 0) part of the composite operator propagator (see our comments after 
Eq. (gOJ)). 

Equating factors in front of ln(— p 2 ) in both sides of Eq. (}4*2"|) . (JSj), arLC ^ 
putting k = 1, 2, . . . m, we obtain the following set of m equations: 



m i r 

Y.Ct\B{l + lik + l) = ^ Jdxx m (l-x)[(l 



xf-'-x'- 1 } 



1=1 

1 

~ 2 



B(m + l,k+l) - 



(m + k){m + k + 1) 



fc= 1,2, ...m.(44) 



8 It corresponds to a theory of free quarks (which interact only via electromagnetic 
forces) . 

9 Here and in what follows we omit terms which are sub-leading in the limit p 2 — * 0. 
The superscript means that a corresponding quantity is calculated in n-th order in 
strong coupling constant a s . 
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Using formulae (|A.2|) - (jA.5|) from Appendix A, we find the solution of these 
equations in the form 10 



£(0) 



c 



(o) 



m 
l- 



l = 1, 2 m 



(45) 
(46) 



where (^) denotes a binomial coefficient. Let us stress, the quantities C m> i 
PB^) satisfy equation for all integer k ^ 1, although for our purpose it 
was enough to take m values of /c. 
One can check that 



m ' 1 (/ + !)(/ + 2) 







(47) 



i=i 



(see Eq. (JA3J)). Note that the RHS of Eq. (|28|) is identically equal to zero 
for k = 1. Indeed, (J J O l,1 )(x,Q 2 /p 2 ) = due to the Furry theorem. On 
the other hand, Eq. (g2D means that (O^O" 1 '')^ ~ [(I + + 2)] -1 . We 
conclude, it is the relation (|4Tj) due to which Eq. is satisfied in the 
leading order in a s . As we will see below, a similar relation is valid for 
as well. 

Now let us consider the next order in a s . The diagrams describing prop- 
agator of the composite operator in this order are presented in Figs. EaJ [50 
and El Using Feynman rules for the composite operators shown in Fig. |HJ we 
obtain the following contribution of the diagrams in Fig. |p"al 



{O n ' h O m ' l )^\p 2 /^ 
2 



-1 



■ n-l 



X 



8^ 

k-l 



C F ln 



—p i 



In 



—p z 



k(k + l) 
B(l + l,k + l) + (k^±l) 



B{l + l,k + l) + J2 B ( l + 1. 3 + 1) 

3=1 



(48) 



10 Let us stress, we didn't demand from the very beginning that C m ^ m should be equal to 
zero for odd m. It is a consequence of the fact that electromagnetic interactions conserve 
P-parity. Remember that DIS structure function F2(x,Q 2 ) is an even function of variable 
x, and its non-zero moments are Fa(n, Q 2 ) — C n ,n (Q 2 /^ 2 ){p | O n ' n |p)(/i 2 ), with n = 2k. 
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where Cp = (N 2 — 1)/2N C , with N c being a number of colors. For the 
diagrams in Fig. I5bl we get the expression: 



(O n ' k O m ' l ) { ^{p 2 /fi 2 ) = (-1) 



47T 3 V ~P' 



In 



~P 



(49) 



r?7 



Both diagrams in Fig. |H1 are equal to zero since they are proportional to 
Thus, the sum of diagrams which make a contribution to the propagator of 
the composite operator is given by 



(O n ' k O m ' l Y 1] (p 2 7 y) = (-lr-^cvin 



in-" 



x 



-p^ 



In 



k(k + l) 



B(l + l,k + l) 



k-l 



+ 2 yU L_ 

^ \k- j k + 1 

3=1 V J 



£(/ + l,j + l) + (/e^/) 



(50) 



Note, that (O n - fe O m '') (1) = for fc = 1 or I = 1. Equation flH| can be 
presented in the form: 



(O n > k O m > l ) {1 \p 2 /fi 2 ) = ^C F \n 

4:71 



X 



(k 
-4^ - - 1 + 
i=2 ^ 



k-l 



+ 2 yp L_ 

^ \k - j k+1 

.7 = 1 V J 



fc(Jfe + l) 

(o n ' j o m 'Y \p 2 /ii 2 ) + (kT±i) 



(51) 



For our further purposes, it is useful to present (|5T)|) in a form which has 
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no explicit symmetry in I and k: 



(O n ' k O m ' l ) {1 \p 2 /fi 2 ) = {-l) n - l p^C F \n 



An 3 



x 



r,2 



— p z 



111 



— p z 



k(k + l) 



B(l + l,k + l) 



k-i , 

2 Er 



— ' V k — j k 



B(l + + 



(52) 



In deriving Eq. ()52|) from Eq. (jBHj) . the formulae from the Appendix A were 
used. Taking into account (pftj)). we find 



m / -2 \ 

J2 C% <0^0"*>«(pV) = (-l^e^^ln In 



x [i + (-iH 



fc-l X 

+ 2 £ ^ 



i=2 



/ 1 + + 1) 



B(m+l,k + l) 



B(m+l,j + 1) 



^ \ k- j k+1 

i=i v J 

+B-i) , (7_" 1 1 ) f(-*Ej-i+ 

z=i v 7 L \ j=2 J 



k(k + l) 



B(l + l,k + l) 



k-l 



+ 2 v J L_ 

~\k — j k + 1 



B(l + l,j + l) 



(53) 



As one can see from (J53J), ££1 Cffi (0"' fc O m '') (1) = for jfc = 1 and m = 1. 
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The expression ()53|) can be simplified and presented as 11 

m / -2 \ 

£ (O^O^Vtf/f) = (-ir^i^m (A^J In 
4 E7 



— 



J'=2 



, 2 2 

1 + + - 

j km 



B(m+l,k + l) 



(m + fc)(m + A; + 1) 



fe-i 



i=2 



k — 3 



B(rn+l,j + 1) - 



+ j)(m + j + 1) 



+ 



2(m - l)(jfe - 1) 
m(m + l)(/c + l)(m + k) 



(54) 



With the use of Eq. (|A.15|) from Appendix A, one can rewrite (|54jl in the 
form which has an explicit symmetry in k, m: 



z=i 
x J 2 



x 



+ 



^ 1 1 , 2 2 

-> -- V - + 1 + - + - 

J'=2 J j= 

B(m + l,k + l) 



k m 
1 



m + k)(m + k + 1) 
V ; B(m + l,j + 1) + J2 : + 1,3 + 1) 

m+fc— 1 

£ 7 + 



fc-l 



(m + k) 



(m + k)(m + k + 1) ^-^ j A; (A; + l)m(m + 1) 



(55) 



Now let us calculate the Green function which contains one composite 



ii 



We used the relation B(m + 1,2) - [(to + l)(m + 2)]" 1 = 0. 
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operator and two electromagnetic currents. The diagrams in Fig. [7a] give 



disc (p+9)2 (J J 0"' fc )« (x, Q 2 /p 2 , p 2 ^ 2 ) 



— 



In 



— 



x 



k{k + l) 



k-l 



x(l - x) k - (1 - x)x k 



k(k + l) 



[a;(l — a;)- 7 — (1 — x)x J 

j=2 



(56) 



The equation (J55j) can be rewritten as follows: 



0:.s 



disc (p+(?)2 (J J O n ' ) ( a '(x,Q /p ,p / fi ) = ^C F \n 



x 



fc(Jfe + 1) 



fc-i 



&s C[p+q? (JJO n > k )^{x,Q 2 /p 2 ) 



k(k + l) 



^disc (p+g)2 (J JO n ^°\x,Q 2 /p 2 )\, 



(57) 



with zero order expression given by 

The contribution from the sum of the diagrams in Fig. [7E] is equal to 



disc (p+g)2 (J J O n > k )^ (x, Q 2 /p 2 , p 2 /fx 2 ) 



_p2 



In 



-p z 



x 



i=2 



XIX 



k-l 



3=2 



k — j k 



x(l — x) 3 — (1 — x)x 3 



(58) 



Note again that Eq. (J58)l can be rewritten in terms of an absorptive part of 
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the 3-point Green function calculated in the leading order in a s : 



disc (p+9)2 (J J O n ' k }^\x,Q' z /p z ,p 2 /fi 2 ) = ^C F \n 



a. 



2tt 



—p z 



2 J2 - disc ( P+9 ) 2 ( J J ° n,fc ) (0) Q 2 /p 2 ) 

3=2 3 ) 

+ E (^3- - 1) disc (^) 2 (^^o^y°)(x,g 2 /p 2 )| 



(59) 



The diagram in Fig.|S]is sub- leading at p 2 — > 0, since it has no singularities 
at p M = 0. The diagrams in Figs. and IHE1 describe a renormalization of 
the electromagnetic current J em which is conserved and, therefore, has no 
anomalous dimensions. As a result, these diagrams are suppressed by the 
inverse of ln(/z 2 /(— p 2 )) with respect to (JoT)j) and (p)%|) . Thus, we have to sum 
only expressions (JBTIj) and 

l 

L y dEP^-ldisC^a (J J O n,k ) ( - 1 \x,Q 2 /p 2 ,p 2 /fJ 2 ) 



,n-l „2 



O.s 



e 9 Q~~3 ^ 111 



_p2 



In 



~P 



x 



i=2 



./ 1 + /- 



£?(m + 1, A; + 1) 



1 



(m + A;)(m + k + 1) 



fc-i 1 



fc+ 1 



B(m + l,k + l) 



1 



(m + j)(m + j + 1) 
1 



k + 



k-l 
J=2 



B(m + 1, j + 1) 



(m + k)(m + k + 1) 
1 



(m + j)(m + j + 1)_ 



(60) 



Making use of equation (jA.7|) . one can easily calculate the sum in the last 
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line of equation (|6U|) and rewrite it in the form: 



i- J dxx m ~ 1 disc (p+9)2 (J J O n ' k } {1) (x,Q 2 /p 2 J p 2 //j 2 ) 



Al 2 2 

-4V--1 + - + - 
^— ' i km 



j=2 



fc-1 1 



fc - j 



5(m+ 1, j + 1) 



P / V ~P 

B(m + l,k + l) - 



(m + k)(m + k + 1) 



i=2 

2(m-l)(A;-l) 
m(m + l)(k + l)(m + fc) 



(m + j)(m + j + 1) 



(61) 



The set of equations for the OPE CP's C^\{Q 2 / fi 2 ) {I = 1, 2, . . . m) looks 



like: 



(0) 



1=1 



1=1 



1 

= 2^/ rfxa;m " ldisc (^) 2 (JJO n ' k )^(x,Q 2 /p 2 ,p 2 /fx 2 ) 


From (J62j) and (|54|) . (|61|) one, therefore, can get m equations, 

E O B(l + l,k + l)=e 2 q ^C F \n (%) 

1=1 m W ' 

k 

- 4 E7 



(62) 



i=2 



2 2 

+ + - 

j km 



B(m + l,k + l) 



(m + k)(m + + 1) 



fe-i 



B(m+l,j + 1) 



1 

(m + j)(m + j + 1) 



2(m-l)(Jfe-l) 



m(m + l)(k + l)(m + fc) 



(63) 
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for m different values of k = 1, 2 ... m. In particular, we obtain for k = 1: 



E 



i 



Z=l 



(/ + !)(/ + 2) 



0. 



(64) 



This equation is similar to formula (|47j) which was obtained in the leading 
order. 

However, in order to find C^\, it is much more convenient to use equiv- 
alent form of the expression in the RHS of Eq. (J63|) For this purpose, one 
should compare with (JHHJ), and then exploit the above demonstrated 
symmetry in k, l{k,m) 

E O B(k + 1, 1 + 1) = C F In g!) 

x { [! + (-!)•»] (-4E--1 + 



i=2 



m(m + 1 
J + l) 



B(k + l,m+l) 



( 1 1 \ 

+ 2 E :--T # + 1. 

z — ' \m — i m + 1 J 

j=i 7 

-eW(T-"i) [(-^r 1 

i=l v 7 L \ J=2 J 



/(/ + !) 



(65) 



Thus, we get m equations (corresponding to k = 1, 2 ... m) for m coefficient 
functions. The solution of these equations (J55J) is rather easy to find, the 
result is 



Cm,m(Q 2 / A* 2 



x[i + (-i; 



i=2 



j m(m + 1) 



(66) 
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and 



2 a s n i 
e ^ CFl11 



f-< 2 



X 



+ 



fc=z+l v 7 v 



m—1 
l-l 

1 

m+1 

m — 1 
fc- 1 



' 1 



1 + 



i=2 



Z(Z + 1) 



k-l 




(67) 



for / = 1, 2, 



1 . 12 We used formulae for summation in I ( (|A.5|) , (|A.6J) ) 
and j ( (jA.lljl ) from Appendix A. In particular, we get from ()67j) (m ^ 2): 



2 «s ^ , 



// 2 



m 1 



1 1 

2 m — 



m + 1 



(68) 



As one can see from (jfifij) . "major" coefficient function C, 
well-known anomalous dimension of O m ' m 



(i) 

m,in 



ryNS 

I III 



^C F 
2tt 



m 1 

4 Et 



J'=2 



m(m + 1) 



is defined by 



(69) 



in spite of the fact that O m,m mixes with all operators O m,i (I — 1, 2. . . . m) 
under the renormalization. We have reproduced the standard expression for 
the coefficient function C m , m (Q 2 //i 2 ), and, simultaneously, have calculated 
"gradient" OPE coefficient functions C m j(Q 2 / fi 2 ) (I — 1, 2. . . .m — 1) in zero 
and first order in a s (see Eqs. ()46|) and (JH7J)). 

Let us emphasize that our operator definition of the coefficient functions 
C m '\Q 2 /fi 2 ) (pTTj) results in a homogeneous renormalization group equation 
for C m ' 1 (Q 2 1 /z 2 ) with respect to /i, in spite of the fact that composite operator 
Green functions (O n)k O m ' 1 ) need an additional renormalization. The details 
are discussed in Appendix B. 

12 It easy to check that our CF's (|S5|> . ljfi7|) satisfy Eqs. not only for k = 1, 2 . . . to, 
but also for arbitrary integer k ^ 1. 
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5 Finite renormalization of composite 
operators and rescaling of coefficient 
functions 



With accounting for equations (pfo"J) . and (jo^|) . (joTj) . the expressions for 
the OPE CF's can be rewritten in the following form: 




for I = 1, 2, . . . , m — 1. 

Let us consider a sum of products of renormalized composite operators 
and corresponding CF's which enter the OPE of two electromagnetic cur- 
rents (j3J. According to (jl2j) . the renormalized composite operator 0^' z (/^) 
is changed under rescaling \i\ — > fi 2 as follows 

0£'(a*?) = E^/^)°^)> ( 72 ) 

where /ij denotes a renormalization scale of the composite operators in the 
MS-scheme, 13 and Z is a matrix of a finite renormalization. Let us emphasize 
that due to Eq. (|B~9j) 

/ k 

(o m > l o n > k )M) = E zM/^l) E ziiAltA) (o m > l 'o n > k ')M)i (73) 

l'=l k'=l 

In the renormalization scheme with subtractions, — /i 2 is off-shell renormalization point 
of Green functions. 
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since additive coefficients in Eq. (|B.9|) do not depent on the renormalization 
scale (jl. 

In the first order of strong interaction, the matrix looks like 



1 + C F In 

47T 



t4 
pi 



1 2 



Zw = I 



n a ° i 

Z7T 



A*2 



i=2 

1 



i(i + i)y 



nij u -v i + i 
o, 



v = i, 

v < i 
v > i 



In a general case we find: 

m ml 

1=1 1=1 l'=l 

m 
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i=i 



where 



k=l 

Thus, we obtain from (J75J) 

Cm,m (^2/^1) = ^m,m Z m (/^//^l); 



(74) 



(75) 



(76) 



(77) 



ctmjQ4/fi) = c m ,i + Yl ^ Zi^l/ti), (78) 

fc=2+l 

for I — 1, 2, . . . , m— 1. It was taken into account that = for k < I. Thus, 
if one changes the renormalization scale /i, the "major" coefficient function 
C m ^ m is simply multiplied by the factor Z mm , while the other CF's mix with 
each other. It is a consequence of the fact that the matrix Z^ has a triangle 
form with zero elements above its diagonal. 

The coefficient functions (7 m ,z (^I/a*?) are normalized so that C m j{l) = 
C m j, where numbers C m j are "bare" coefficient functions which were ob- 
tained for the case with no strong interactions (see Eqs. (|4"5|). (jUJ). We 
conclude from explicit expression of Z in the first order in a s (J74j) that our 
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perturbative result, equations (f7U|) and (|7Tj). is a particular case of general 
formulae (J77D and (|7%|) . 14 



6 Conclusions 

As was shown above, the OPE coefficient functions can be expressed in terms 
of the Green functions of the corresponding composite operators without ex- 
plicit use of the elementary (quark and gluon) field operators. We would like 
to specially stress that our general expression (p?T|) for the OPE coefficient 
functions holds in any renormalization scheme in contrast with previous pre- 
scriptions (see, for instance, jlOj). We believe that such a representation 
is useful in proving general properties of the coefficient functions and for 
non-perturbative calculations (e.g. in lattice QCD). 

We would like to stress that if one assumes the light-cone expansion in 
the framework of axiomatic approach (see, for instance, Ref. ^H] ) ? then our 
formulae do not lean on perturbative expansions at all and deal with v.e.v. 
of T-products of local Heisenberg operators. 

In this paper we limited ourselves by the QCD nonsinglet operator only. 
The singlet case needs a special treatment and will be considered in a forth- 
coming paper. 



Appendix A 

In this Appendix we have collected formulae which are needed for calculating 
sums in I and j presented in the text and getting compact expressions. Let 
us first consider summation in index I. For m ^ 0, we have the relation 



D 

1=0 



m 



l + k + 1 
l + l 



1 



/ + 1 m + k+1 



(A.l) 



The quantity (?) is a binomial coefficient. 15 In the following equations, m ^ 1 
is assumed: 



■m-l / \ 

£m>"; 

i=i v 7 



:-ir+ 



r(2-m) 



(A.2) 



14 After an obvious replacement mI/Mi ~ > Q 2 /^ 2 - 

15 The beta-function is given by the equation [B(n + 1, k + = (n + k + 1) { a \ k ) 
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m—1 



m — 1\ 1 1 



E(-i)\7- 1 "jT = -^ |1 + ( - iri - 

1=1 v 7 



(A.3) 
(A.4) 



l)(/ + 2) (m + l)(m + 2) 



[l + (-in, (A.5) 



-2 y ( - J— - -J— 

j' \m — t m+l 



(/ + !)(/ + 2) 



(m + l)(m + 2) 



m 1 



i=2 



m(m + 1) 



(A.6) 



m— 1 



, /m - 1\ (I + A; + 1 



m+l 



i-i J V i + i 

m + fc + lV 1 1 



m + l J m+l (m + k){m + k + 1) 



, (A.7) 



- f_ir +1 / m + ^ + 1 ^ 'I + k L_ c a 8) 

1 m + l y m (A; + l)(m + £; + !) m + k A ' ; 



gv^x^rr 1 



Z(Z + 1) 

= / _ [ m + k + lN \ 1 * ca 9) 

1 m + l y m(m + l) (fc + l)(m + fc + 1) 
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Now let us consider summation in index j (k,l ^ 1 is assumed every- 
where) : 



fm + k — j — 1\ 1 (m + k — 1 
pi V m ~ 1 ) 3 



k 



m+k—l m—1 

. i=i J i=i J 



(A.10) 



TO— 1 Z-l 1 / 1 \ -1 / • i \ i 

El \ 1 (m — 1\ \ m — 7 — 1\ 1 
7-E- = ( / _ 1 £' 



=1 J 3=1 J V 



(A.ll) 



•7 = 1 " J = l - ' i = l X / J 

Note that the last equation is a particular case of Eq. (jA.lOJ) for 1 < I < m. 
Let us present other useful formulae: 



^Xfk + j + v - 1 



3=1 



J 



fc + Z + l /fc + Z + l 

k V Z 



1 



(A.12) 



X 



fc + Z + 1 



' 1 ^ (k + A 1 



.i^i J 



3=1 



3/3 



(A.13) 



E 

3=1 



fc + j 



- 1 



fc-1 1 
3=1 J 



- 1 



(A. 14) 



1 



/7 i ' i 1\ "I 1 



3=0 



3 



k + l + 1 



l + k + 1 



k-j 



k 



3 =2 J . 
-1 k ' 



3=2 



.(A.15) 



The last relation is a consequence of Eqs. ()A.13j) and ()A.14J) . It demonstrates 
that the sum in (|A.15J) is symmetric under replacement I ^ k. 
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The last useful formula contains sums in / and j: 

m—l / -,\ /i , i , t\ — l 1 k 



m - 1 \ I + k + 1 



l-l 



l + l 



1 



1 + 



J / J 



m + + r 

m+1 J m + 1 



'fc-i 

E 



m + j — 1^ fc — j + 1 

j 



(k-j)(m + j) 



+ (-irE 



m + j\ 1 
J / J 



(A.16) 
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Appendix B 



The goal of this Appendix is to derive a renormalization group equation (see, 
for instance, [15]) for the CF's stating from our formula (pTT|). Let us rewrite 
it in the form: 16 

m 1 } 

C m ,i (00)% = — / dxx m ~ l disc (p+g)2 (JJO) k R (B.l) 

l=i m { 

(see our brief notations in the end of Section |2J). The renormalized quantity 
in the RHS of Eq. (jB.lj) is given by 

Y,(J\n)(n\JO n > k ) R) (B.2) 

n 

where Yl n means a sum m a complete set of elementary (quark and gluon) 
states \n). Let us underline that each state |n) should include at least one 
quark- ant iquark pair (in non-singlet state). It is well-known that a matrix 
element with an insertion of one composite operator is multiplicatively renor- 
malized (see, for instance, [14J). Therefore, one has the following relation: 

(J\n) R = Zj(J\n) v . (B.3) 

Since electromagnetic current is conserved, its renormalization constant, Zj, 
does not depend on the renormalization scale, [djd\i)Zj = 0. , 

Now we will show that the matrix element {n\JO n,k ) r in (jB.2|) is also mul- 
tiplicatively renormalized. Consider arbitrary diagram G which contributes 
to a matrix element with an insertion of two composite operators, J and O k . 
We assume that divergences of all sub-diagrams of G are already removed. 
Let G to have (p + 2) external lines (p 0), and q internal vertexes. The 
index of this diagram, u(G), is defined by [H] 

q 1 p+2 

u(G) = < aX + 4 - ^ + 2 ) + ("J - 4 ) + - 4 )' ( B - 4 ) 

i I 

Here id™ "* is a maximal index of internal vertex of a type i, and ri is a power 
of a polynomial corresponding to an external field of a type I (r q = 1 for 

16 Hcrc and in what follows the subscript ft/m means that a corresponding quantity is a 
renormalized (unrenormalized) one. 
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a quark line). The quantities uj is the dimension of the electromagnetic 
current, while UlQn,k IS the dimensions of our composite operators. Note that 
u max q f or a u QCJ3 vertexes, and Uj — 3. 

As for u ™,k, it is formally equal to (& + 2) (only derivatives with respect 
to quark line momenta have to be taken into account, not total derivatives). 
In particular, u; n,i = 3. Let us show, however, that effectively u n,k does 
not depend on k, and Let {ri} (i = 1, . . .p + 4) to be a set 

of external momenta (p,q including), while {ki} (i = l,...m) to be a set 
of loop momenta of G (with m being a number of loops). After using the 
Feynman parametrization and linear transformations from {ki} to {/»}, an 
analytical expression for the diagram can be written in the form: 

„ k-l 

1(G) = H dx j 5(l-J2^s)J2 a ^ 

j J s r=0 

x f[[-FTTnM r Z^y'K > ( B -5) 
t£J ( 2?r ) D(x l ,llr i r ) 

where the Feynman parameters. The change of variables {k{\ — ► {/»} 

is done so that the denominator in Eq. ()B.5jl depends on U only via If. 
Notice, the momentum l m is a linear combination of r^, and k m , ingoing 
quark momentum for the vertex O n,k . The numerator in Eq. (jB.5|) is a sum 
of polynomials of (hlj), (hrj) and {riTj). Each polynomial includes at least 
one scalar product of the vector n, namely, (Un) or (t^ti). 
After integrating in d D li, . . . d D l m _i, we get 



» fc-i 

i(G) = n d Xj 5(i-j2^)j2 b 



[X.; 



r=0 



(2ir) DK D(x u P,r irj ) 

(we denoted l m = I). Note that dim/) — dimP = 5. Let us analyze possible 
terms in P which is a polynomial of I 2 , (Irj), (r^rj), and depends linearly on 
(In) or (fjn). The term (/n)(/ 2 )^ dimP_1 ^ 2 results in zero (for even r) or an 
expression proportional to = (for odd r) in the dimensional regular- 
ization. 17 All other terms proportional to (In) drop for the same reason (if 

17 Similar arguments are used under calculations of the anomalous dimensions of the 
composite operators (see, for instance, Ref. |S]). 
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r ^ 1) or give finite integrals which converge in the ultra-violet region. One 
can verify that terms proportional to (r^n) give non-zero convergent integrals 
only for r = 0. 

Thus, we have to put io Qn ,k = cJo n l = 3 in (jB.4|) . In its turn, it means 
that irf{G) ^ —1, with uj(G) = —1 for p — 1. In other words, there is no need 
to introduce a new counterterm : J(x)0(y) : 5(x — y), and, consequently, 

k 

(n\JO n ' k ) R = Zj x J2( Z ^' WO n ' k ')u. (B.7) 

k'=l 

It follows from (IB~3i (IB~7l) that 

k 

disc (p+9)2 (JJO) k R = Zf ^(Z- 1 )^, disc (p+g)2 (JJO%. (B.8) 

k'=l 

This result is in agreement with perturbative QCD calculations from Sec- 
tion m 

Now let us turn to the propagators of the composite operators. The 
renormalization properties of Green functions with insertions of more than 
one composite operator were studied in details in Ref. [T7]. In particular, 

/ k 

(00)% = ^(Z- 1 )!, ^(Z- 1 )*, [(00)$' + ff] . (B.9) 
u=i k'=i 

The divergent coefficients fg are, in general, non-zero even in the free theory, 
as it is the case for our quark operators Q n,k . 

The renormalization matrix of the composite operator Z depends on a 
renormalization scale /i, while the coefficients f B do not. The renormalized 
Green function (00)% is also /i-dependent, but it does not depend on the 
regularization scale /L All physical (measurable) quantities should be, of 
course, /2 and //-independent. 

As one can see, Eq. (jB.9|) which relates renormalized and unrenormalized 
propagators has an additive term. Nevertheless, renormalization group equa- 
tions for the propagators and for CF's have no additive terms. 18 Indeed, act- 
ing on both sides of Eq. (jB.l|) by the operator \idjd\i = fid/dfi + /3(g s )d/dg s , 

18 However, the renormalization group equations are not homogeneous if some of the 
composite operators has non-zero vev [17] . In our case, (O n ' k ) = 0, and O n ' k do not mix 
with the identity operator under the renormalization. 



30 



and taking into account Eqs. (|B.8J) . (|B.9|) . we obtain: 



1=1 



d ~ 



E d ™,l (00)% - C mtl J24 (00)%' 



l'=l 



where 



7 



z- 1 



(B.10) 



(B.ll) 



is the matrix of anomalous dimensions of the composite operators. Equa- 
tion (|B~T0|) can be represented in the form: 



m , m 

i=i afi u=i 



0. 



(B.12) 



Since (|B.12j) is valid for arbitrary integer k ^ 1, we derive a set of renor- 
malization group equations for the coefficient functions C m j (I — 1, 2, . . . m): 



d 

i^-r- C m ,i — 2_s ^ ^ m ''' 



0. 



(B.13) 



In particular, the coefficient function C mim , which is relevant for DIS, obeys 
the following closed equation: 



j _ m \ / > 



0. 



(B.14) 



Remember that 7 is the matrix of the anomalous dimensions of non- 
singlet composite operators which enter light-cone OPE (jHJ): 



d 



dfj, ' 



(B.15) 



i'=i 



This renormalization group equation demonstrates again that composite op- 
erators mix under the renormalization. After a diagonalization of the renor- 
malization matrix Z (|T3*|) . one obtains new operators O m ' 1 = Jj v =i A\< O m ' 1 ' 
which are multiplicatively renormalized. The matrix A has a triangle form, 
with detA 7^ 0. The anomalous dimensions of the operators O m ' 1 are deter- 
mined by diagonal elements of Z. The one-loop values of 7; = 7/ are given 
by Eq. (jHSD- 
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The analysis of renormalization properties of the composite operator 
propagator is simplified, if we choose a light-cone axial gauge n^A^ = 
in which only diagrams shown in Fig. El contribute to (00) kl . Correspond- 
ingly, the renormalization of the composite operator in this gauge is given by 
diagram in Fig. ^2 The blob in these figures is a sum of all possible QCD 
diagrams (with a disconnected part included). A specific character of the 
diagrams in Fig. HH is the following: they can be divided in two parts by 
cutting two quark lines. It enables us to derive the following relation: 

(00)S = A$ v [e, {tf/p 2 Y] Z% l \e, fj V \ (B.16) 

where the matrix Zb(£, p?/p 2 ) is a bare quark loop with an insertions of two 
composite operators: 

Z B l (e,jl 2 /p 2 ) ~ - e (J^y B(k + 1,1 + 1) (B.17) 

(terms which are non-singular at p 2 —>■ are omitted in (jB . 1 7|) ) . Thus, the 
coefficients fjg in ()B.9|) . which are needed to regularized the Green functions 
of two composite operators, look like 

fS~~A$ v {e,l)BW + l,l' + l). (B.18) 

The relation (jB.16|) is confirmed by our calculations in perturbation theory 
(see, for instance, one-loop (jUj) and two-expressions (JBTjl for (00) r). 
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Figure 1: The propagator of the composite operator, (<f>n k 0m i) ■> an d the 
Green function ((f) 2 2 k ) in free scalar field theory. 
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Figure 2: Feynman rules for the composite operators O™'^ ^ n in the leading 
(zero) order in strong coupling a s . 
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Figure 3: Feynman rules for the composite operators 0^ ,k ^ n in the first order 
in strong coupling a s . 
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Figure 5b: The diagrams for the propag ator (O n > k O m ' l )W (continued). 




Figure 6: The diagrams which give zero contribution to the propagator 

lnn,kQm,l\ (1) _ 
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Figure 7b: The diagrams for the matrix element (JJO n ' k )^ (continued). 
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Figure 8: The diagrams which give non-leading contribution to the matrix 
element (JJO n ' k )^ at p 2 -> 0. 
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Figure 9a: The diagrams which give contribution to the renormalization of 
the electromagnetic current J em . 
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Figure 10: The diagrams which give contribution to the propagator of the 
composite operator in the axial gauge. 
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Figure 11: The diagram which describes the renormalization of the composite 
operator in the axial gauge. 
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